Spin-wave theory for dimerized ferromagnetic chains 
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Abstract. We describe a Peierls dimerization which occurs in ferromagnetic spin chains at 
finite temperature, within the modified spin-wave theory. Usual spin-wave theory is modified 
by introducing a Lagrange multiplier which enforces a nonmagnetic state at finite temperature. 
It is shown that this method gives results in excellent agreement with the density-matrix 
renormalization group applied to transfer matrices for dimerized ferromagnetic chains. We study 
bond correlation functions (S^ • St+i) and explore the characteristic features of dimerization in 
the specific heat. 



1. Introduction 

It is a well established fact that a coupling between electronic and lattice degrees of freedom can 
drive structural instabilities of electronic systems. This effect is induced by a gain in electronic 
energy surpassing the cost in elastic energy due to the lattice distortion. For instance, in the 
Peierls instability, a static lattice distortion turns a one-dimensional (ID) free electron system 
into a band insulator pQ. A similar effect has been observed for antiferromagnetic (AFM) spin 
chains coupled to phonons where a gain in magnetic energy drives the so-called spin-Peierls 
transition [2j. A further example for a Peierls instability, where a ferromagnetic (FM) spin- 
chain shows a periodic modulation (dimerization) of the magnetic exchange in a certain finite 
temperature region was observed in YVO3 [3] and has recently been discussed theoretically [HE]. 
This effect is induced by a coupling of the spin and orbital degrees of freedom [6] . 

In this paper we highlight the differences of the spin-Peierls transition in AFM and FM 
spin chains and discuss the dimerized FM spin chain within the framework of modified spin- 
wave theory (MSWT). We show that this method yields results in excellent agreement with 
numerical data [5]. In particular, we consider the specific heat and bond correlation functions. 
In addition we give an analytical expression for the free energy at low temperatures and small 
dimerization allowing for a quantitative understanding of the occurrence of dimerization in FM 
spin chains. 



2. Spin-Peierls effect in spin chains 

First, we investigate under which circumstances a Peierls dimerization may occur in FM spin 
chains. It will be also instructive to discuss briefly the Peierls dimerization in AFM spin chains 
to highlight the differences between the FM and AFM case. 

Our starting point is a coupling of the spin chain to lattice degrees of freedom. If the Peierls 
gap exceeds the phonon frequency it is permissible to neglect the phonon fluctuations. Within 
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Figure 1. Phase diagrams obtained from equation (1) using TMRG for: (a) the dimerized 
AFM S = 1/2 Heisenberg chain — the inset shows S(T) with K/J = 2, (b) the dimerized FM 
S = 1/2 Heisenberg chain — the inset shows 5(T) with K/\ J\ = 0.1. 



this approximation the Hamiltonian reads 

N 

H = J E i 1 + c- 1 )**} s * • Vi + ™? , (i) 

where we have defined S = 2gu/(Jao) and K = KJ 2 aQ/(Ag 2 ). Here J is the exchange constant, 
Sj a spin S operator at site j, N the number of sites, and K the effective elastic constant 
with E e \ = NKu/2. The dimensionless parameter 5 is given by the exchange constant J, the 
spin-phonon coupling constant g, the atomic displacement u, and the lattice constant ao- We 
note that in spite of the fact that the model ([T|) is ID, the static, mean-field treatment of the 
three-dimensional phonons allows for a finite temperature phase transition if 5{T) is treated as 
a thermodynamical degree of freedom. 

For the ID Heisenberg antiferromagnet it was shown that -Emag ~ —5 A ^ [7j. Hence the cost 
in elastic energy due to a distortion of the lattice is outweighed by the contribution stemming 
from the magnetic part of the Hamiltonian. This leads to a spin-Peierls transition for every value 
of K at sufficiently low temperature. The phase diagram obtained within the density-matrix 
renormalization group applied to transfer matrices (TMRG) [8] is shown in figure (TJ^a) for the 
dimerized AFM 5 = 1/2 Heisenberg chain (pQ). The phase transition is of second order. The 
inset shows the parameter 5 as a function of temperature. 

In figure HJb) the phase diagram for the dimerized FM chain obtained by TMRG is shown. 
We observe that a dimerized phase exists only at T > and only if K/\ J\ < K C /\J\ ~ 0.118. We 
also report a tricritical point (TCP) at (T TC p/\J\, K TC p/\J\) (0.696,0.116). For T < T TC p 
and T > Ttcp the transition will be first and second order, respectively. This is shown in the 
inset of figuredjb) where we can see that at the upper phase boundary 5(T) evolves continuously. 

3. Modified spin-wave theory for dimerized chains 

For the case of a FM chain, we apply a Holstein-Primakoff transformation onto the magnetic 
part of Hamiltonian (pQ). Since the unit cell is doubled for the dimerized chain we have to 
distinguish between magnons on two sublattices. Hence we write 

S i = \I 2S ~ a )<* a 3* a j« ' S Jcc = a \ a \/ 2S ~ a ) a a ja . S l* = S ~ a L a i-« > ( 2 ) 



where a^ a and aj Q are bosonic annihilation and creation operators at site j. The index a S {e, o} 
refers to the sublattice for j being an even (e) or an odd (o) index, respectively. Retaining only 
the lowest order terms we obtain a bilinear Hamiltonian which can be diagonalized by means of 
a Bogoliubov transformation. The diagonalized Hamiltonian reads 



#mag = H + J2(oJ k a\a k + u+@l/3 k ) , (3) 
k 



with H = -\J\S 2 N and the magnon branches uo k = 2|J|S'(1 =t \J cos 2 k + 5 2 sin 2 k) which are 
shown in figure [2(a). 

The idea of Takahashi's MSWT [§] is to retain Hamiltonian © while enforcing the Mermin- 
Wagner theorem. This is achieved by minimizing the free energy subject to the constraint that 
the magnetization at finite temperature vanishes. To this end a Lagrange multiplier is introduced 
serving as a chemical potential. Applying MSWT for the dimerized FM chain we obtain the 
constraint 

S 4S{" b K) + «bK + )} , (4) 

k 

with n^{u k ) = {exp[(u;^ : — n)/T] — l} -1 being the Bose factors. The free energy per site in this 
approximation reads 

/ = -| J\ S + - i ln t 1 + nB K~)] + ^ [1 + n B K + )] } • (5) 

k 

For the uniform case (5 = 0) it reduces to the free energy obtained by Takahashi for the ID 
Heisenberg ferromagnet [9]. 

The free energy can be obtained analytically within the MSWT framework in the limit 
ts = T/{\J\S(1 - 5 2 )} < I. To this end we expand equation (JH) [10J to obtain 



tg ^1 + ct\fT$ + ^a 2 t$ + . . . \ 



(6) 



with v = —4fiS 2 /T, and a = )/5y / 7r. Inserting this result into equation ([5]) we obtain 

/ = -\J\S- T{C(|)^ /2 /(2v^) + 0(t s )}. Hence the gain in magnetic energy is proportional 6 
as is the elastic energy, explaining the existence of a critical value K C /\J\. 

Due to the alternating interaction between the spins, nearest-neighbor correlation functions 
(NNCF) will also exhibit alternation in strength. Hence Ass = {&2j • S27+1) — {S%j ■ S2J-1) 
serves as an order parameter for the dimerized chain. In calculating the corresponding NNCF 
within MSWT it is essential to also retain quartic bosonic terms which leads to 

(S 2j • S 2i±1 ) = (l 53 [«bK) - n B K + )] ft(S)\ , (7) 



with f k (5) = (cos 2 k ± 5 sin 2 fe)/( v cos 2 fc + 5 2 sin 2 fc) . In figure [2(b) results from MSWT for 
A55 as a function of 5 at various temperatures are compared to those obtained within TMRG. 
The agreement is good up to T/\J\ ~ 1 and even the right limit is captured by MSWT for the 
fully dimerized case. However, corrections 5 = 1 — e are predicted to be of order e 2 for e<l 
whereas perturbation theory exhibits corrections of order e. 

With equation ([7]) it is easy to determine the internal energy per site reading 



u = \ {(1 - 6) (Szj-i ■ S 2J ) + (1 + S) (S 2j ■ Sjy+i>} 



(8) 




Figure 2. MSWT results for S = 1: (a) Magnon dispersions lo^ (lower branches) and lj"^ (upper 
branches) for various values of 5 £ [0,0.9]. (b) Ass as a function of 5 for various T. The 
lines (symbols) denote the MSWT (TMRG) results, (c) Specific heat per site as a function of 
temperature for various 5. 



Using this result the specific heat per site reads c = c + + c with 

c± = £ (-cr)n B [1+71B (Wfc/)] (w& - /* + T^) /± {n B ( w -) - n B K + )} /± (9) 

k.k', 

ae{±} 

In figure [2jc) we show the specific heat per site as a function of temperature for various 
dimerizations. For 5 = we find one broad hump which still is retained if 5 <C 1. However, 
at a higher dimerization we observe that the two magnon branches give distinct contributions 
leading to a sharp low temperature peak and a hump at higher temperatures stemming from 
well separated magnon excitations, see figure E^a). The excitations of the lower magnon band 
u7 may occur at low T whereas those of the higher band uj~£ contribute only when T is high. 

4. Conclusions 

We have discussed the differences in the occurrence of a spin-Peierls transition in AFM and FM 
spin chains. For the FM chain we have applied MSWT and found that the magnetic energy gain 
~ —T 3 / 2 5 2 directly competes with the cost in elastic energy. The dimerization order parameter 
obtained from MSWT is found to be in excellent agreement with TMRG data. Finally, we have 
shown that dimerization leads to a qualitatively different behaviour of the specific heat. 
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